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Abstract 



A dataset of 126,501 spiral galaxies taken from Sloan 
Digital Sky Survey was used to analyze the large- 
scale galaxy handedness in different regions of the 
local universe. The analysis was automated by using 
a transformation of the galaxy images to their ra- 
dial intensity plots, which allows automatic analysis 
of the galaxy spin and can therefore be used to an- 
alyze a large galaxy dataset. The results show that 
the local universe (z<0.3) is not isotropic in terms 
of galaxy spin, with probability P <~ 5.8 • 10~ 6 of 
such asymmetry to occur by chance. The handed- 
ness asymmetries exhibit an approximate cosine de- 
pendence, and the most likely dipole axis was found 
at RA=132°, DEC=32° with la error range of 107° 
to 179° for the RA. The probability of such axis to 
occur by chance is P < 1.95 • 10~ 5 . The amplitude 
of the handedness asymmetry reported in this paper 
is generally in agreement with Longo, but the statis- 
tical significance is improved by a factor of 40, and 
the direction of the axis disagrees somewhat. 

Keywords: Cosmic parity, Galaxies, Cosmolog- 
ical Principle, Computational Astrophysics, Dipole 
Axis 



1 Introduction 

Assuming that the universe is isotropic, it is ex- 
pected that in a sufficiently large sector of the uni- 
verse the number of galaxies that rotate clockwise 
will be roughly equal to the number of galaxies that 
rotate counterclockwise. However, recent evidence 
suggests that the local universe does not follow that 
expected balance, and show that the ratio between 
clockwise and counterclockwise galaxies in some re- 
gions is significantly different than 1:1, introducing 
galaxy handedness asymmetry. 

Longo [l| analyzed the sense of rotation of 15,158 
spiral galaxies imaged by Sloan Digital Sky Survey 
(SDSS) and showed cosmic parity violation and a pos- 
sible cosmic dipole axis. The analysis was based on a 
relatively small dataset of galaxies classified manually 
by four undergraduate students, so that the number 
of galaxies in each 30° Right Ascension (RA) slice 
ranged between to 3512 galaxies. Redshift of the 
galaxies was smaller than 0.085. The experiment was 
done using a web-based user interface that mirrored 
half of the galaxy images, so that the analysis was 
not biased by certain preferences of the human read- 
ers who classified them. A dipole axis was detected 
with RA=217°, DEC=32° with a probability of oc- 
curring by chance of 7.9- 10~ 4 [1]. 

Asymmetry of clockwise and counterclockwise 
galaxies was also observed in the Galaxy Zoo dataset 
[21 ■ However, since the classification was done man- 
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ually by a static graphical user interface, the asym- 
metry could be attributed to certain preferences of 
the Galaxy Zoo participants who used the web-based 
system. 

The sense of rotation of a galaxy is a gross metrics 
that is relatively easy to measure, and it is not af- 
fected by atmospheric effects or hardware inaccuracy. 
However, analyzing large datasets of spiral galaxies 
manually requires significant labor, and might be bi- 
ased by human errors or reader preferences. In this 
study a computer analysis was used to determine the 
handedness of spiral SDSS galaxies, and the results 
were analyzed in the light of the Cosmological Prin- 
ciple and the cosmological isotropy assumption. 

2 Image dataset 

The data used in the experiment are galaxy images 
collected by SDSS [3[, and classified by Galaxy Zoo 
4]. The first full Galaxy Zoo data release contained 
891,552 galaxies with redshift, classified manually by 
Galaxy Zoo participants. The galaxies that were 
classified manually as spiral galaxies were selected, 
and the Ganalyzer method described in Section [3] 
was used to reject the edge-on galaxies to provide 
a dataset that contains only spiral galaxies that their 
handedness can be determined. That provided a 
dataset of 352,705 spiral galaxies with identifiable 
sense of rotation. It should be noted that the galaxy 
population is not uniformly distributed in the night 
sky covered by SDSS, and some RA ranges have little 
or no galaxies. 

The galaxy dataset of Galaxy Zoo was determined 
based on the human perception, and using a static 
user interface. Therefore, while the results of Galaxy 
Zoo provide a "super-clean" galaxy dataset, they 
might be biased by human behavior and preferences. 
To avoid possible human bias, another dataset that 
was used in the experiment was the SpecObj view of 
SDSS DR7, which is based on the SpecObjAll table, 
but excludes bad and duplicate data. The galaxies 
in the SpecObj view were classified automatically us- 
ing the Ganalyzer method [B|, providing a dataset of 
345,599 spiral galaxies that are not edge-on. Declina- 
tion of the galaxies was between ~-11.2° to ~70.3°. 



The Galaxy Zoo dataset is available on-line through 
the Galaxy Zoo Data Release 0], and the SpecObj 
view is available through SDSS CAS server. 

3 Image analysis method 

The handedness of each spiral galaxy was determined 
by transformation of the images to their radial inten- 
sity plots, as done by the Ganalyzer method 0, 
To transform the galaxy image to its radial inten- 
sity plot, first the foreground pixels of the galaxy 
are separated from the background using the Otsu 
global threshold Q- Then, the center of the galaxy 
is detected, and the radius is determined by the max- 
imal distance between the object center and any fore- 
ground pixel 

The radial intensity plot is a 360x35 image, such 
that the value of the pixel {x, y) is the median value 
of the 5x5 windows around the pixel at image coor- 
dinates {O x + sin{6) ■ r, O y — cos{9) ■ r) in the galaxy 
image, where 9 is the polar angle (in degrees) and ris 
a radial distance. Intuitively, the radial intensity plot 
is an image of the radial intensities at different dis- 
tances from the galaxy center. Each horizontal line 
in the radial intensity plot is then smoothed using a 
median filter with a span of 50 pixels. 

Figure [T] shows the original galaxy image and a 
transformation of the radial intensity plot such that 
the Y axis is the intensity and the X axis is the polar 
angle. As the figure shows, in an image of a spiral 
galaxy the peaks are expected to shift due to the spi- 
rality of the arms, while in galaxies that are not spiral 
the peaks are expected to align in a near-straight line. 
The radial intensity plot transformation of galaxy im- 
ages is described more thoroughly in 

Once the radial intensity plot is generated, the 
peaks of each curve displayed in FigureQ]are detected 
using a peak detection algorithm §, and the slope 
of the line created by the peaks can indicate whether 
the galaxy rotates clockwise or counterclockwise [Ij]. 
The slope is measured by first finding the peaks in the 
radial intensity curves displayed in Figure [TJ and sep- 
arating the peaks into groups such that each group 
contains peaks of the same arm of the spiral galaxy 
Figure [2] shows the peaks detected in the radial 
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Figure 1: Galaxy images (left) and the transforma- 
tion of the radial intensity plots such that the Y axis 
is the intensity and the X axis is the polar angle 

intensity plots of Figure Q] 



Figure 2: The peaks detected in the radial intensity 
plots of the elliptical galaxies of Figure [1] 

For each group of peaks, the peaks are ordered by 
their Y coordinate in the radial intensity plot. Then, 
each peak at X coordinate x in the radial intensity 
plot is compared to the X coordinate of its next peak. 
If the X coordinate of the next peak is smaller than 
x then the counter I is incremented. If the X coordi- 
nate of the next peak is greater than x, the counter 
r is incremented. If after searching through all peaks 
I > 3r, it is determined that the galaxy rotates clock- 
wise. If r > 31 then the galaxy rotates counterclock- 
wise. All other cases are considered undetermined 
handedness, and are excluded from the rest of the 
analysis. 

Using a set of 120 spiral galaxies taken from 
Galaxy Zoo and used in 0, and available for 
download at 



http://vfacstaff.ltu.edu/lshamir/downloads/ganalyzer 
the method was tested for its accuracy and error rate. 
All galaxies had clear spirality, and the handedness, 
of each galaxy was classified manually by the author. 
The results showed that when using the criteria 



described above of I > 3r and r > 31 to determine the 
handedness, the handedness of 42.5% of the galaxies 
was in agreement with the manual classification 
performed by the author, while the remaining 57.5% 
of the galaxies had undetermined handedness. When 
using I > 2r and r > 21 as criteria the percentage 
of galaxies with determined handedness increased 
to ^68.3%, but with the sacrifice of ^6.7% error 
rate among those galaxies that their handedness was 
determined by the algorithm. When using I > r 
and r > I as criteria the determined handedness 
increased to 96%, but the error rate also increased 
to 15%. 

Once the handedness of the galaxies is determined, 
the handedness asymmetry of a certain population of 
galaxies can be determined by A = fej, where A 
is the asymmetry value, L is the number of galaxies 
rotating clockwise, and R is the number of galaxies 
rotating counterclockwise. 

The accuracy of the Ganalyzer method is sensi- 
tive to redshift and magnitude when compared to 
human classification. However, Galaxy Zoo and the 
SpecObj view are not random samples, and contain 
bright objects that can be processed by Ganalyzer 
as described and analyzed in [1]. In lower redshift 
of <0.07 Ganayzer was able to determine the hand- 
edness (I > 3r or r > 31) of ~42% of the galaxies 
in the SpecObj view that were classified as spirals, 
while in the redshift ranges of 0.07-0.14 and 0.14- 
0.21 the detection rate was ~38% and ~36%, respec- 
tively. While the rate of the detection of the hand- 
edness is different at different redshift ranges, it is 
expected that in each redshift range the detection 
rate of clockwise rotating galaxies will be equal to 
counterclockwise rotating galaxies. The dependence 
on redshift can also be compared to manual classifica- 
tion of galaxies. While at redshift of 0.05 Ganalyzer 
classification into spiral and elliptical galaxies was in 
agreement with the manual classification in ~89% of 
the cases, at redshift 0.1 it was ~87%, and ~82% 
and ~74% at redshifts 0.15 and 0.25, respectively 0]. 
/Tiik^xh4iMa)^3s,.zi|3i be also attributed to the classifi- 
cation accuracy of the Galaxy Zoo citizen scientists, 
which drops as the redshift gets higher The red- 
shift dependence of Ganalyzer accuracy 0, [ltj can 
affect the magnitude of the signal of the asymmetry, 
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but any inaccuracy is expected to be equal to both 
clockwise and counterclockwise directions, and there- 
fore the direction of the asymmetry should not be af- 
fected. The Pearson correlation between the degree 
of spirality measured by Ganalyzer [f| and the ap- 
parent magnitude computed using ~60,000 galaxies 
with redshift <0.3 provided a weak Pearson correla- 
tion of -0.036 [nj, with P< 10~ 8 . That experiment 
showed that in the datasets of Galaxy Zoo and the 
SpecObj view, which contain brighter galaxies com- 
pared to the galaxies of SDSS DR7, there was merely 
a weak dependence between magnitude and the abil- 
ity of Ganalyzer to detect spirality. Galaxies in which 
spirality was not detected were rejected and therefore 
did not have an effect on the analysis. 

To test the response of Ganalyzer to changes in 
the brightness of the galaxies the experiment was re- 
peated such that the I > 3r and r > 31 criteria were 
used, and the intensity of the foreground pixels of 
the galaxies was scaled down by different factors. As 
Table Q] shows, the percentage of galaxies with deter- 
mined handedness remains stable until the intensities 
of the foreground pixels are scaled down by more than 
0.25. However, even when Ganalyzer has difficulties 
in determining the handedness, it does not missclas- 
sify left handed and right handed galaxies. The sig- 
nal of the radial intensity plots weakens as the galaxy 
gets dimmer compared to its background, and there- 
fore less peaks are detected when the galaxy is dim- 
mer. However, the slope of the peaks determined by 
Ganalyzer [5j cannot be inverted because less peaks 
are found. 
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Table 1: Classification accuracy when the intensity of 
the foreground pixels of the galaxies is scaled down 



4 Results 

The method described in Section[3]was applied to the 
two galaxy datasets described in Section [5] The anal- 
ysis was based on galaxies that met the handedness 
criteria described in Section [3l so that Galaxy Zoo 
and the SpecObj view provided 126,501 and 121,414 
spiral galaxies with determined directionality, respec- 
tively. 

Figure |3] shows the handedness asymmetry of the 
SpecObj galaxies defined in different RA ranges and 
three redshift ranges, which are <0.085 (y, 0.085- 
0.17, and 0.17-0.3. The standard error was computed 
based on normal distribution using 1/y/N, where N 
is the total number of galaxies in the region and red- 
shift range. The RA ranges [60,90] and [270,300] are 
empty due to the very low number of galaxies (70 and 
0, respectively) in these RA slices. 
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Figure 3: Galaxy handedness asymmetry in different 
RA ranges and redshift ranges using SpecObj galax- 
ies. The dotted line shows the cosine curve multiplied 
by -0.04 and starts at 132°, to show cosine depen- 
dence as will be discussed later in the paper 

As the figure shows, different RA and redshift 
ranges show different galaxy handedness asymme- 
tries, suggesting that the local universe is not 
isotropic. Partial correlation can be observed be- 
tween the redshift ranges in some of the RA re- 
gions, but the high standard errors in the different 
RA/redshift bins do not allow statistically signifi- 
cant conclusions about such correlation. Also, as de- 
scribed in Section [3] there is dependence between the 
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rcdshift of the galaxies and the ability of Ganalyzer 
to detect their spirality and determine their handed- 
ness. 

Figure |4] shows the galaxies handedness asymmetry 
defined in different regions using both Galaxy Zoo 
dataset and the galaxies from SpecObj. 
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Figure 4: Galaxy handedness asymmetry in different 
RA ranges 

The differences between Galaxy Zoo and SpecObj 
can be explained by the perceptional bias of the 
Galaxy Zoo dataset, in which the human-based clas- 
sifications are biased towards left-handed spirals by 
~10%. However, it is clear that the two datasets 
provide the same profile of galaxy handedness asym- 
metry. The strongest asymmetry of the SpecObj 
data was 0.028±0.00683, found in the RA range 
[120°, 150°]. Assuming that the probability of a 
galaxy in that RA range to rotate clockwise is equal 
to its probability to rotate counterclockwise, given 
the standard error the probability to observe asym- 
metry of 0.028 by chance is ~ 4.3 • 10~ 5 . Considering 
just the four RA slices with more than 20,000 galax- 
ies in them (the fifth most populated RA slice has 
just 9907), the probability that one or more of them 
have asymmetry equal or greater than 0.028 and one 
or more other slice have asymmetry of 0.16 (as ob- 
served in 180°-210°) or greater given the standard 
errors in these RA slices is^5.8-10~ 6 . 

As a test for bias, the analyses were repeated with 
the dataset described in Section [2] such that ran- 
domly selected half of the galaxies were mirrored. 
The handedness asymmetry values in the different 



RA ranges in the SpecObj and Galaxy Zoo data were 
smaller compared to the original data, and in all RA 
ranges the asymmetry was within the standard er- 
ror, indicating that no systematic error was present. 
The largest handedness asymmetry measured with 
the mirrored data was 0.009±0. 01683. 

A real spiral spin preference would present itself 
as a dipole with a cos </> dependence, where is the 
space angle between the position of the galaxy and 
the preferred axis [l[. The RA, DEC and asymmetry 
values of each galaxy were used to find the most likely 
dipole axis using cosine dependence [l[ . Each galaxy 
i in the dataset was assigned with {—1,1} handedness 
di, and the <fi angle was computed for each galaxy as 
the difference between the angle of the galaxy and 
the angle of the candidate dipole axis. The d ■ | cos 4>\ 
was then fitted to cos0 for all possible RA, DEC 
combinations in an exhaustive grid search such that 
the declination range was [-10°, 70°], within the sky 
covered by SDSS, to find the most likely dipole axis 
with which the galaxy handedness asymmetries have 
cosine dependence. 

To deduce the most likely dipole axis, random 
handedness values were assigned to the galaxies, and 
the random handedness were fitted to cos 4> for each 
right ascension and declination axis, and repeated 
1000 times for each axis to compute its mean and 
standard deviation of the x 2 computed with random 
galaxy handedness. Then, the x 2 computed for the 
same axis with the observed galaxy handedness de- 
termined by the Ganalyzer method as described in 
Section [3] was compared with the mean and a of the 
X 2 computed for that axis using the 1000 runs of ran- 
dom galaxy handedness values. The axis that had the 
highest difference between the random and observed 
galaxy handedness in terms of a was found around 
RA=132°, and DEC=32°. Figure |5] shows the dif- 
ference in terms of a between the observed \ 2 an d 
the mean % 2 of 1000 attempt of fitting randomized 
handedness values to dipole axes with varying RA 
and declination 32°. 

As the figure shows, the most likely dipole axis is 
at RA=132°, and the la error range for the dipole 
is between RA 107° and 179°. With - 4.27ct the 
probability to have such axis by chance is P < 1.95 • 
10~ 5 , showing that it is likely that galaxy handedness 
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tion range, respectively. The asymmetry values at 
different declination ranges and RA slices are dis- 
played in Figure |U 
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Figure 5: a between the mean % 2 computed with 
1000 attempts of fitting all galaxies with randomized 
handedness values to dipole axes with varying RA 
and declination 32° 



asymmetry form a dipole axis in the local universe 
(z < 0.3). Handedness asymmetry in the 30° RA 
range around the most likely dipole [117°, 147°] is ~- 
0.017, and the asymmetry is ^-0.027 in the same RA 
range where z<0.085. The declination of the most 
likely dipole axis is in agreement with the most likely 
dipole axis proposed by [l[, but the RA of 132° is 
different from the RA=217° [l[. However ,the error 
range of the most likely dipole was ±35° [H, within 
the lrj statistical error range reported here. 

As the graph also shows, cosine dependence was 
observed in low (z<0.085) and higher (0.085<z<0.3) 
redshift ranges, and the dipole axis is consistent in 
both ranges. The signal for the dipole is far stronger 
for z<0.085. The weaker signal at the higher red- 
shifts can be attributed to the ability of Ganalyzer 
to accurately detect spirality, which decreases as the 
redshift of the galaxies gets higher [5|, ll(| . The most 
likely dipole axis in z<0.085 is at RA=136°, with 
probability of ~ 7.19 • 10~ 5 of having such cosine de- 
pendency by chance. 

The handedness asymmetry was also analyzed for 
different declination ranges. The declination ranges 
that were used are [-10°, 10°], [10°,30°], and [30°,50°] 
of the SpecObj view galaxy dataset, which have 
31,943, 46,482, and 33,006 galaxies in each declina- 
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Figure 6: Galaxy handedness asymmetry in differ- 
ent RA and DEC ranges. DEC ranges are [-10°, 10°], 
[10°,30°], [30°,50°] 

The limited declination ranges covered by SDSS 
do not allow comparison of declination ranges in op- 
posite hemispheres, and not all declination and RA 
ranges of SDSS have galaxies in them. The highest 
number of galaxies in these declination ranges are in 
the RA range of [150°,210°], which contained 69,986 
galaxies. As Figure [6] shows, the galaxy handedness 
asymmetry in that RA range increases with the dec- 
lination, while the asymmetry in other RA slices such 
as [210°, 240°] the galaxies in the higher declination 
are biased toward counterclockwise rotation. 

5 Conclusion 

The results of the experiment show that galaxy hand- 
edness asymmetry is different in different RA ranges 
within z<0.3, and the asymmetries exhibit a cosmo- 
logical dipole axis. According to these observations, 
the local universe is not isotropic, meaning that the 
observed physical characteristics of the universe are 
different in different directions of observation. These 
results are in agreement with [l|. The most likely 
dipole axis was found at RA=132°, DEC=32°. This 
axis is somewhat different from the most likely dipole 
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axis reported by at RA=217°, DEC=32° with an 
uncertainty of ^35°, but the difference is within the 
error ranges. 

The limited range of declination in the SDSS 
dataset as well as the absence of data in certain right 
ascension ranges did not allow comparing all opposite 
hemispheres for handedness asymmetry for the pur- 
pose of fully profiling the characteristics of the hand- 
edness asymmetry in the local universe. It should be 
noted that the experiment is limited to the local uni- 
verse, although the range of z<0.3 exceeds far beyond 
the scale of a galaxy supercluster. Higher-resolution 
analysis of the cosmic parity violation in the more 
distant universe will be possible when the Large Syn- 
optic Sky Survey (LSST) starts operating, providing 
a much larger galaxy dataset and seeing deeper space. 

The experiment is based solely on data acquired by 
Sloan Digital Sky Survey, and under the assumption 
that these raw data are not biased. A possible source 
of bias might be the automatic galaxy detection al- 
gorithm that is part of the SDSS pipeline, that might 
have a preference to galaxies of some certain direc- 
tionality. This, however, conflicts with the observa- 
tion that some of the opposite hemispheres have op- 
posite asymmetries. In any case, even in the case that 
such bias exists, it is expected to be consistent for 
all SDSS data, so that the asymmetry values might 
be shifted, but the asymmetry profile itself should 
not change. Another possible weakness is that galax- 
ies from different RA were taken at different times, 
and any change made to SDSS or its processing al- 
gorithms during that time could potentially lead to 
a biased dataset, and different galaxy asymmetries 
at different redshifts. It should also be noted that 
the SpecObj view and the Galaxy Zoo dataset con- 
tain bright objects, and are not random samples of 
the galaxies in DR7. A bias in the selection of these 
objects can also affect the results. 

Clearly, this study does not offer explanations nei- 
ther to the nature of asymmetry at the cosmologi- 
cal scale, nor to the observation that the asymmetry 
changes at different regions. 
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